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d\ ' 

O ' Abstract 

o ■ 

I We consider a system with a single base station communicating with multiple users over orthogonal channels 

while being assisted by multiple relays. Several recent works have suggested that, in such a scenario, selection, 
i.e., a single relay helping the source, is the best relaying option in terms of the resulting complexity and overhead. 
However, in a multiuser setting, optimal relay assignment is a combinatorial problem. In this paper, we formulate 
a related convex optimization problem that provides an extremely tight upper bound on performance and show that 
selection is, almost always, inherent in the solution. We also provide a heuristic to find a close-to-optimal relay 
assignment and power allocation across users supported by a single relay. Simulation results using realistic channel 
models demonstrate the efficacy of the proposed schemes, but also raise the question as to whether the gains from 
relaying are worth the additional costs. 
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I. Introduction 



^ . In distributed wireless systems wherein each node possesses only a single antenna, relays can be used 

^ to provide spatial diversity and combat the impact of fading. Relaying has been an extremely active 

^ ' research area, especially since Sendonaris et al., in [1], proposed the idea of user cooperation wherein 
O 

Q\ mobile users cooperate by relaying each others' data. Many cooperation schemes have now been studied. 



e.g., [l]-[4]. The work in [2] and [3] proposed repetition-based cooperation schemes including fixed 



^ amplify-and-forward (AF) and decode-and- forward (DF) using orthogonal channels (time/frequency slots). 
' In networks with multiple relays, the traditional strategy has been to let all the relays forward their 
messages to the destination. However, having relays transmit on orthogonal bands is bandwidth inefficient. 
A proposed alternative is to use distributed space-time codes (DSTC) [3]; however, this requires symbol 
level synchronization, which is difficult to implement over a distributed network. It has recently been 
shown that most of the benefits of cooperative diversity can be achieved with minimum overhead if a 
single 'best' relay cooperates with the source. This scheme is referred to as selection cooperation [5], [6] 
and has now been investigated in various contexts [5]-[9]. 

In the case of a single source-destination pair, choosing the best relay is fairly straightforward and 
solved for both DF [5], [6] and AF [7] relaying. In both cases, the best relay is the one that contributes 
the most to the output signal-to-noise ratio (SNR). The selection gets significantly more complicated in 
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the more practical case of multiple information flows [6]. Because a relay must now divide its available 
power between all flows it supports, a relay that is best for a single flow may not remain the best overall 
and relay selection becomes a combinatorial problem. In [6], the authors present ad hoc approaches to 
approximate the optimal solution with limited complexity, without addressing resource allocation. 

In relay networks, an independent research theme is that of resource allocation, including power 
allocation, e.g., [10], [11] amongst many. Optimal allocation makes best use of the limited available power 
resources. However, to our knowledge there has not been any work jointly considering relay assignment 
and power allocation in the context of multiple source, multi-relay, especially cellular, networks. 

Our system model comprises a single base station communicating to multiple users being assisted by 
a few dedicated relays. The users are to be assigned to the relays. The relays have limited power which 
must be divided among the users they support. Relaying in the context of a cellular wireless network 
has received limited attention [12]. In Section |Inl we develop an optimization problem for optimal relay 
assignment and power allocation at the relays. We try to answer the question, what relay assignment and 
power allocation scheme maximizes the sum rate and what scheme maximizes the minimum rate to all 
the users! Obtaining solutions to these requires exponential complexity. 

The main theoretical contribution of this paper is in Section Hill where we derive upper bounds to the 
rates and show that these bounds form a convex optimization problem for both figures of merit. We use 
the resulting Karush-Kuhn-Tucker (KKT) conditions to illustrate why the bound is tight and then derive a 
simplified, tight, lower bound. In Section |Wl we simulate a cellular network, using the COST-231 model, 
to study the performance gains in a relay assisted network over a traditional single base station system. 
Interestingly, while the gains are significant, the results leave open the question of whether these gains 
adequately compensate for the additional infrastructure costs of a relay-assisted cellular system. 

In terms of the available literature, our formulation is similar to tone assignment in orthogonal frequency 
division multiple access (OFDMA) systems. This is the problem of assigning users to tones to maximize 
a certain metric, subject to the constraint that no tone is assigned to two different users. In [13], Wong 
et al. solve the problem of minimizing the total transmission power in a multi-user OFDM system. The 
frequency spectrum is divided into discrete frequency bins and a convex relaxation technique is used to 
solve the discrete bin-assignment problem. Rhee and Cioffi [14] solve the problem of assigning users to 
tones on the downlink of a OFDMA based communication system with the objective of maximizing the 
minimum rate to users. 
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While our approach is similar to that taken in these papers; the roles of 'user' and 'tone' in these papers 
is played, respectively, by 'relay' and 'user' here. However, there are some significant differences, most 
importantly in the power allocation step. There a single power allocation is required across all tones, here 
each relay must meet a power constraint Furthermore, our problem formulation allows us to analyze the 
conditions under which the bounds are tight unlike the other works, wherein the authors state that the 
bound gets tight as the number of tones approaches infinity, but prove this for only the two user case. 

This paper is organized as follows. In Section |Ill we describe the system model in some detail. In 
Section [nil we then formulate the optimization problem and the upper bound to each of the two rates 
and illustrate why the bounds are tight. In Section |IVl we illustrate this through simulations and use the 
bound to analyze the performance of relay assisted cellular networks. The paper wraps up with some 
conclusions in Section IVl 

II. System Model 

Our system model consists of a cellular network with a single BS, communicating with K users, and 
assisted by J relays, as shown in Figure dJ Each of the users is assigned an orthogonal channel, over 
which the BS-to-user and the relay-to-user communications take place. The users are frequency division 
multiplexed, although the results here also apply to the case of time division multiplexing. The relays in 
the system are fixed wireless terminals, installed solely to aid the BS-user communication. The relays use 
the DF protocol with the same codebook as the transmitter. 

The communication between the BS and a user happens over two time slots. In the first time slot the 
BS transmits, while the relays and the user try to decode the message. In the second time slot, one of 
the relays, chosen a priori, re-encodes and then transmits the information it has decoded in the first time 
slot. The user uses the messages received in the two time slots to decode the transmitted information. 

Suppose that user k (denoted as d^) is allotted to relay-j (r^). For a system as described above, the 
maximum rate at which the BS can communicate with the receiver with the help of the relay is [6] 

Jrf^ = mm{lsr^,Isrjdk) , (1) 
Isr, = ^ logs (1 + SNR,|/l,,^ P) , (2) 
Isr.d, = ^log2 (1 + SNR,|/l,rfj2 + SNR,ajfc|/l,^.rfj2) , (3) 

where, SNR^ and SNR^ are, respectively, the ratios of the transmit power at the BS (denoted as s) and 
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relay to the noise power at the receiver, hgr^ is the channel between the BS and relay j, denoted by rj, 
similarly K-d^ is the channel between relay rj and destination d^. Finally, ajk is the fraction of the total 
relay power used to communicate with user k. The factor of 1/2 accounts for the fact that the BS-user 
communication happens over two time slots. Isr^ is the rate at which the source can communicate with 
relay -j while Isrjdk is the maximum rate at which the source can communicate to user k with the help 
of the relay. Equation ([T]) ensures that both the relay and the user can decode the message. 

The channels between the BS, relays and users are modeled using the COST-231 model as recommended 
by the IEEE 802. 16j working group [15]. The model includes the path loss, large-scale fading (a log-normal 
variable) and small-scale fading modeled as Rician random variable for line-of-sight (LoS) communication 
and Rayleigh random variable for non-LoS communication. When the BS and relays are both placed at 
some height above the ground, the fading has a LoS component. The existence of this component is 
crucial since it suggests that all relays will be able to decode a source codeword; hence the factor limiting 
the overall rate is the second term of Eqn. ([T]), Isr^d^^ significantly simplifying the problem at hand. 

III. Problem formulation and solution 

As described in the previous section, every user is assigned one of the J relays. This paper deals 
with optimizing this assignment to maximize two metrics of interest, the sum rate to all the users and 
the minimum of all the rates. In maximizing the sum rate (equivalently the average rate), the objective 
function is 

K K 

^^d^ = Y^ "^i^ {^sr{d^). Isrid,)d,) , (4) 
k=l k=l 

while in maximizing the minimum rate, the objective function is given by 

mm{/dj = mmjmin J^^(rf^)rfJ} , k = l,...,K, (5) 

where, in both cases, r(dk) is the relay assigned to user k. 

In practice, the number of users, K, will be much larger than the number of relays, J. Hence, a single 
relay will likely be required to support multiple users, and to meet its power constraint it must divide 
its power amongst these users. Thus, our objective is now two fold, one, finding the relay assignment 
scheme, and two, once the assignment is done, distributing powers at each of the relays amongst the users 
it supports. 

To formulate a tractable problem, in this paper we investigate simplified versions of the above problems. 



As mentioned earlier, in a cellular network, the data rate bottleneck is the compound source-relay- 
destination channel, the second term in Eqn. ([T]). We assume that 

hr, > Isr.dk Vj, k, (6) 

and hence min {lsr{d^)Jsr(d^)d^) = Lr{d^)d,- 

In Section |Wl we justify the validity of this assumption. Note that in spite of the assumption, the 
solution is not immediate. The fact that the relays divide their power amongst the users they support, 
makes the relay assignment an integer programming problem with the attendant exponential complexity. 

A. Max sum rate 

The sum rate measures the maximum throughput delivered by the base station. For the sake of brevity, 
let Cfc represent SNR^I/is^^p and pjk represent SNR,.|/irjrfj.p, j = 1,2, . . . , J. Let ajk be the fraction of 
the power of relay-j, used to communicate to user k. The optimization problem maximizes the sum rate 
to all the users subject to two constraints: only a single relay helps each user and each relay must meet 
a power constraint. The formal optimization problem is, therefore. 



X^2^°^M 1 + + J^Pifcttj-fc J , (7) 
fe=i V j=i / 

such that Wk, ajkaik = 0, j I, j,l e {1,2, . . . , J} , (8) 



max R = max 



K 



5^«,fc = l V J, (9) 



k=l 



ajk > 0, (10) 

where the objective function assumes the relay uses the same codebook as the source. Equation ([8]) enforces 
the selection rule allowing only one ajk term to be non-zero for all relays. The remaining two constraints 
force the power allocated to be positive and to meet a power constraint. The constraint in Eqn. ^ can 
also be written as an inequality constraint, J2k=i ^jk ^ 1)^J- The solution to the optimization problem in 
either case is the same because the objective function is an increasing function of the powers, ajkS. We 
cannot use the usual gradient based methods to maximize the objective function in Eqn. (|7]). Note that an 
inherent assumption is that the BS has knowledge of the parameters that define the problem. How this 
information is conveyed to the BS is beyond the scope of this paper. 

The solution to the optimization problem in Eqns. (r7])-(fT0l) is complicated by the constraint in Eqn. ([8]). 
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An exhaustive search to find the solution would involve the following: for a given relay assignment, 
solving J water-filling problems corresponding to the power allocation at each of the relays. We need 
do this for every relay assignment and find the maximum of them. Each of the users can be assigned to 
any of the relays, hence, all possible relay assignments must be tested. Doing so is impossible for 
realistic values of J and K. We therefore explore tractable approximate formulations. 

The objective function of the optimization problem in Eqns. (r7])-(fT0l) is concave and the constraints, 
other than the one in ([8]), are affine. Our strategy to solve the optimization problem in hand is to ignore 
the constraints in Eqn. ([8]) and maximize the objective function subject to the power constraints alone: 

min - XI 9 logs 1 + Cfc + (11) 

{",fc} k=l ^ \ i=l / 

K 

such that X«jA: — 1=0 V j, (12) 

fc=i 

-oijk < 0. (13) 

Since we ignore a constraint, the solution so obtained will be an upper bound to the maximum sum rate 
achieved by selection. Note that since this minimization problem is now convex, solving this simplified 
problem is fairly straightforward, e.g., using interior point methods [16]. The computational complexity 
involved in solving the optimization problem is polynomial in K and J, and the problem is, hence, 
tractable for practical values of K and J. We use CVX, a package for specifying and solving convex 
programs [17], [18]. 

We now proceed to show that although we did not impose the selection rule explicitly, the solution to 
the optimization problem has the property that, for most k, aj^aik = 0, j 7^ l,j, / G {1, 2, . . . , J} . This 
means, when the power is optimally allocated, most users receive power from only one of the relays. 
Tightness of Bound: The objective and the constraint functions are differentiable and the constraint 
conditions satisfy Slater's condition [16]. To show this, consider one possible choice for the power 
vectors, ajk = 1/K. This meets the sum power constraint with equality and the constraint on positive 
power with strict inequality. For a convex optimization problem with differentiable objective and constraint 
functions, which also satisfy Slater's condition, the solution to the optimization problem satisfies the KKT 
conditions [16]. 

Let us characterize the set of solutions to the optimization problem. For the sake of clarity, we start 



with the case with two relays. In such a case, the Lagrangian of the minimization problem is given by 

K K 



>C({aifc,a2fc};{Afc},{Afe},z/i,i/2) = -i? - A^aifc - ^ A^asfc 

fc=l k=l 

^aik-U +V2\y2a2k-U , (14) 

vfe=l / Vfc=l / 

where and \\,k = 1,2, . . . , K are the Lagrange multipliers associated with the constraint on positive 
power, and ui and 1/2 are the Lagrange multipliers associated with the constraint on the total power at the 
two relays. Any solution to the optimization problem satisfies the KKT conditions, which are, 

— + A^ = ^i, \laik = 0, \l>0, (15) 

— + A^ = ^2, Xla2k = 0, \l>0. (16) 

Now suppose for some z G {1, 2, . . . , K}, an and a2i are both non-zero, then the conditions X]aii = 
and Xja2i = dictate that \} and A^ are both zero. From the KKT conditions it follows that 

— = — . (17) 

Pli P2i 

Similarly if aij and a2j are both non-zero for some j G {1, 2, . . . , K}, then 

^ = ^. (18) 

Pi] P2j 

Unless pii/p2i = Pij/p2j, Eqns. (flTl) and (fTSi) cannot simultaneously be true. In the current problem, pjk 
represent the powers of the channels between the relays and the users. If they are independent continuous 
random variables, as is the case with the wireless channels, then the probability that pii/p2i = Pij/p2j 
is zero. Hence, when the power is optimally allocated, at most one of the K (aik,a2k) pairs has two 
non-zero entries and K — 1 of the pairs have at most one non-zero entry. This indicates that the selection 
rule, {aika2k = 0,VA;), which we did not explicitly impose, is true for at least all but one of the K 
users. Hence, the solution obtained by ignoring Eqn. ([8]) comes quite close to the solution to the original 
optimization problem in Eqns. (|71)-(fT0l). 



For the case of three relays, the KKT conditions are 

Pik 

1 + + Y^^i^iPik^ik 

P2k 





^1, 


A^aifc — 


0, 


Al>0, 


(19) 






\2 

'^k'^2k — 


0, 


XI >0, 


(20) 


+ ^l = 


^^3, 


Xla-sk = 


0, 


• Al > 0, 


(21) 



1 + Cfc + 'E^i^iPik^ik 
P3k 

1 + Cfc + Yli=l Pik^ik 

where, for A; G {1, 2, ... , K}, Xl, A| and are the Lagrangian multipliers associated with the constraint 
on positive power, and ui, V2 and z/3 are the Lagrangian multipliers associated with the constraint on total 
power. In the solution to the optimization problem, we wish to find the maximum number of triplets {an, 
Oi2i, Oisi), in which more than one entry is non-zero. We do this by analyzing different possibilities for the 
solution. Suppose that in the solution, for some i, (an, a2i, a^i) are all non-zero (user i receives power 
from all relays), then, 

l>l Vo 

— = — = —. (22) 

Pli P2i PSi 

Now, for some j, if aij and a2j are non-zero, then, along with Eqn. (|22l) . this would imply that Pii/p2i = 
Pij/p2j, which occurs with probability zero. Hence, if the solution to the optimization problem has one 
triplet with all non-zero entries, then all other triplets can have only one non-zero entry, i.e., selection is 
imposed on all other users. 

Now suppose that in the solution, for no i, (an, a2i, «3i) are all non-zero. Without loss of generality, 
suppose for some j, aij and «2j are non-zero, and for some k, a2k and a-^k are non-zero, then, 

Z/i 1/2 V2 ^3 

— = — , — = (23) 

Plj P2j P2k P3k 

These two equations imply that in all other three-tuples {aik, a2k, osa,), only one of the entries is non-zero. 
This is because, if for some /, an and a-^ are non-zero, then, (|23]) would imply, pn/psi = P2jPii/P3jP2i, 
which occurs with probability zero. Hence, for the case of three relays, at most two of triplets can have 
more than one non-zero entry. Like with the case of two relays, when the power is allocated optimally, 
the selection rule is followed in most of the triplets. 

Generalizing this to J relays, when the power is allocated optimally, at most J — 1 of the J-tuples 
(aik, a2k, • • • , Oijk) can have more than one non-zero entry. This indicates that if ^ J — 1, as expected 
in practice, a large fraction of the users are guaranteed to receive power from only one relay. 

To summarize, we have shown that the power allocation matrix, [ajk]jxK is sparse. Most of the rows 



of the matrix have only one non-zero entry. If all the rows of the matrix had at most a single non-zero 
entry, then we would have obtained the solution to the optimization problem given by Eqns. (|7])-(fT0l). A 
simple heuristic to find that solution, then, is to explicitly impose selection: assign users receiving power 
from multiple relays to the relays that allot the maximum power. 

r{dk) = rm if m = aj:gmax{ajkPjk} (24) 

j 

If there are multiple relays which allot the same maximum power, assign the user to any one of them 
arbitrarily. Once this relay assignment is done, J water-filling problems can be solved for the power 
distribution at each of the relays. However, we can also re-use the power allocation vector derived from 
the earlier step. Construct the matrix [a'jf^]jxK as follows: for each /c G {1, 2, . . . K}, 

<^'mk = <^mk a'jk = VjV ^- (25) 

The matrix of the power allocation vectors [a'jf^jxK meet the constraints given by Eqns. ([8]) and (flOl) and 
satisfy J2k=i '^'jk — Ij^J- It hence a lower bound to the solution to the optimization problem given by 
Eqns. (r7])-(fT0l). We avoid a second round of optimization because, as we shall see, the upper and lower 
bounds are already indistinguishable. 

B. Max sum rate with a minimum rate constraint 

Maximizing the sum rate does not ensure any fairness with respect to the distribution of power. In a 
cellular network, a more practical metric might be maximizing the sum rate while guaranteeing a minimum 
rate to each user. Formally, the optimization problem is. 



^l^^og^ll + Ck + ^PjkOijkU (26) 

k=l V 3=1 J 

such that Vfc, ajkaik = 0, j 7^ / G {1, 2, . . . , J} (27) 



max R = min 

{ajk} {ajk} 



K 



J]a,fc-1 = V J, (28) 



k=l 



< 0, (29) 
i?fc - ^ log2 ^1 + Cfc + ^Pjfeajfc^ < 0, (30) 
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where is the rate guaranteed to user k. Suppose we ignore the constraint given in Eqn. (|27]) . the 
Lagrangian of the resulting optimization problem, for the case of J = 2 relays is given by : 

K K 



£({aifc,a2A:};{Afc},{A^},i^i,z/2,{7fc}) = - R - ^ Xlaik - ^ Xla2k 

k=l k=l 

+Z/1 ^ "ifc - 1 + a2fc - 1 



,k=l / \k=l 



(^Rk - ^ log2 ^1 + Cfe + Y^^Pjkajk^ j , (31) 



where A^ and \l, k = 1,2, . . . , K are the Lagrange multipliers associated with the constraint on positive 
power, ui and 1^2 are the Lagrange multipliers associated with the constraint on the total power at the 
two relays, and 7^, k = 1,2, . . . , K are the Lagrange multipliers associated with the constraint on the 
minimum rate. The solution, if it exists, satisfies the KKT conditions, which are, 

^ , , ^''1^ ~2^'^ r + A^ = ^i, A^ai, = 0, A^ > 0, (32) 

2 (1 + Ck + 2^i=iPikO:ik) 

+ Xi = U2, Xia2k = 0, Xi>0, (33) 



P2fc(2-7fe) , ,2 ,, x2„ n \2 



2 (1 + Cfc + Y.i=lPikOiik) 

Ik i^k - ^ log2 ^1 + Cfc + Y^Pjkajk^ ) ^ ^' ^^^^ 

Suppose for some i E {1,2,..., K}, an and a2i are both non-zero, then the conditions Xjau = and 
X^O!2i = dictate that Xj and A^ are both zero, and from the KKT conditions it follows that 

— = — . (35) 

Pit P2i 

Similarly if aij and a2j are both non-zero for some j E {1,2,..., K}, then 

Ui V2 

— = —, (36) 

hence, like the case with the max sum rate metric, when the power is optimally allocated, at most one 
of the K {aik, a2k) pairs can have more than one non-zero entry. To generalize this result, with J relays 
at most J — 1 of the J-tuples (aik, a2k, ■ ■ ■ , ajk) have more than one non-zero entry. The additional 
constraint on minimum rate given in Eqn. (l30l) does not alter this property of the solution. 

A lower bound to the solution of the optimization problem, can be formulated like before. Explicitly 
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impose selection, by assigning users receiving power from multiple relays to the relays that allot the 
maximum power, 



r{dk)=rm if m = SiTgmax{ajkPjk} ■ (37) 

j 

This relay assignment has to be followed with solving J water-filling problems to meet the constraint of 
a minimum rate to each user. Note that it is possible that the simplified optimization problem is feasible 
where as the original optimization problem in Eqns. (|26l) -(l30l) is not. It is also possible that solving the J 
water-filling problems to compute the lower bound might be an infeasible optimization problem. In these 
cases, the bounds are not meaningful. However, these scenarios occur very rarely. 

C. Max-min rate 

We will now consider a third metric, the minimum rate to each user. The optimal power allocation 
ensures that each user receives the same data rate. The optimization problem maximizes the minimum 
rate to all the users subject to two constraints: only a single relay helps each user and each relay must 
meet a power constraint. The optimization problem is. 



max mm , 



^ log2 ^1 + Cfe + ^Pjfc«jfc^ I (38) 

such that V/c, ajkOtik = 0, j ^ I, i,l e {1,2, . . . , J} (39) 

K 

$^a,fc = l V J, (40) 

k=l 

a.-fc >0 V k,j. (41) 



K 



As before, ther than Eqn. (1391 ), the optimization problem in Eqns. (|381)-(|4TI) is concave: the objective 
function is concave and the remaining constraints are affine [16]. As before, we ignore the constraints 
given in Eqn. (l39l) and maximize the objective function subject to the power constraints alone: 



max R = max min 



^log2 ^l + Cfc + ^Pjfeajfcjj, (42) 

K 

such that ajk = 1 V j, (43) 

k=l 

ajk>0 V k,j. (44) 
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Note that the objective function given by Eqn. (|42l) is not differentiable. To analyze this problem, we 
formulate an equivalent optimization problem with differentiable objective and constraint functions. 

The logarithm function is a monotonically increasing function of its argument, and hence, maximizing 
the minimum of logarithm functions is same as maximizing the minimum of the arguments of the logarithm 
function. Also, for any positive real numbers xi, 0:2, ... , Xn, 



mm 

i 



[xi] = ( max 



1 



and max{xj} = lim | 
Using these relations, the objective function can be reformulated as 



X, 



i=l 



mm — 



mm 



i/i- 



where 



1 + Cfc + Y.'j=iPjkOijk- 



Consider the optimization problem for some finite /, 



mm [ — ) 
such that 



mm I > , 

\h (Ri) 

K 

^ a^-fc - 1 = 



i/i 



V J, 



k=l 

-otjk < 



V k,j. 



i/i 



(45) 



(46) 

(47) 
(48) 



Again we show that the solution to the relaxed problem leads to selection in most cases. Note that 
while these conditions are the same as derived earlier for the max-sum rate, the approach to this getting 
here is very different. The optimization problem given by Eqns. (|46l)-(|48]) is also a convex optimization 
problem [16]. To show that when the power is optimally allocated, most users receive power only from 
one of the relays, let us characterize the set of solutions to the optimization problem. For the sake of 
clarity, we again start by looking at the case with J = 2 relays. In such a case, the Lagrangian of the 
minimization problem is given by 

K ^ \ K K 



C{{aik, a2k} , {Afc}, {Afc}, z/i, z/2) 



k=l 



k=l 




(49) 



where Rl = (1 + + pikOtik + P2kOi2k)', ^\ and A^, k = 1,2, ...,K are the Lagrange multipliers 
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associated with the constraint of positive power given by Eqn. (|48l) ; and and z/2 are the Lagrange 
multipliers associated with the total power constraints given by Eqn. (|47] ). The KKT conditions, which 
must be satisfied, are 



K K 



^aifc = l, ^a2fc = l (50) 

fc=l k=l 

-«ifc<0 Vj,A; (51) 

- + z/i = 0, A^aife = 0, A^>0, Vfc, (52) 

-tefc - A^ + z/2 = 0, A2a2fc = 0, A2>0, V/c, (53) 



where 



'''' 

Now suppose for some z G {1, 2, . . . , 7^}, an and Q;2j are both non-zero, then the conditions X]aii = 
and A^Q;2i = dictate that X} and A^ are both zero. ^From the KKT conditions it follows that 

^ = ^. (54) 

Pli P2i 

As discussed earlier, no other pair (ai^, a2k) can have two non-zero entries. Note that this property of the 
solution is true for all /. Therefore, the solution to the optimization problem with the objective function 
given by Eqn. (|45l) . and constraints given by Eqns. (I43l)-(l44l). also has this property. The objective function 
given by Eqn. (l45l) is a reformulation of Eqn. (l42l) . and, in turn, the solution to the optimization problem 
given by equations Eqns. (|42l) - ((44)) also has the aforementioned property for the case of two relays. 
Therefore, for the case of two relays, the solution obtained by ignoring Eqn. (|39| ) comes quite close to 
the solution to the original optimization problem given by Eqns. (|38l) - (|4TI) . 

To generalize this result, the solution to the max-min optimization problem given by Eqns. (I42l)-(l44l). 
for the case of J relays has at most J — 1 of the J-tuples (aik, a2k, • • • , ajfc) with more than one non-zero 
entry. 

The construction of a heuristic to the solution of the optimization problem given by Eqns. (|38])- (|4TI) . 
is, as before: assign each user receiving power from multiple relays to that relay from which it receives 
maximum power. If there are multiple relays which allot the same maximum power, assign the user to 
any one of them arbitrarily. Once this relay assignment is done, if required, J max-min power allocation 
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algorithms are solved for the power distribution at each of the relays. 

D. Independent codebooks at the relays 

In the previous section, relay selection and power allocation was done for the case when the transmitter 
and the relays use the same codebooks to encode the messages. The results can also be extended to 
the case when independent codebooks are employed at the source and the relays. Using independent 
codebooks results in higher rates [3], however, decoding of the source and relay messages is significantly 
more complex compared to the case of repetition coding [19]. When the source and the relays employ 
independent Gaussian codebooks, the optimization problem to maximize the sum rate to all users, similar 
to the ones given in equations (IT])-®, is given by: 



5^|^^°S2(l + Cfc) + ^log2 + |, (55) 

such that Vfc, ajkaik = 0, j 7^ /, j, / e {1, 2, . . . , J} (56) 



max R = max 



K 



^ ajk = 1 V J, (57) 

k=l 

ajk > 0. (58) 

It is not hard to show that other than the constraint given in Eqn. (|56l ). the optimization problem is 
a concave maximization problem, and like before, solving it gives an upper bound to the sum rate. The 
heuristic which also serves as a lower bound can also be constructed from it. A max-min optimization 
problem can also be formulated in a similar manner. 

IV. Numerical Results and Discussion 

In this section we verify the validity of the assumption in Eqn. © and present the results of simulations 
to illustrate the tightness of the bounds developed in the previous section. We compare the performance of 
three cases: the baseline scenario uses a single-input single-output system (SISO) with a single antenna at 
the BS and user and relaying is not used. The alternative is a system with a single antenna at the BS and 
J relays with a single antenna each. The last system considered is a multiple-input single-output (MISO) 
system with J + 1 antennas at the BS and a single antenna at each user. In comparing these cases, all 
other system parameters, e.g., number of users, total power and bandwidth, remain constant. 
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TABLE I 

Parameters used in COST231 model 



Parameter 


Value chosen 


Parameter 


Value chosen 


BS height 


50m 


Rooftop height 


30m 


Relay height 


50m 


User height 


1.5m 


Frequency 


IGHz 


Road orientation 


90 degrees 


Building spacing 


50m 


Street width 


12m 


Transmit power 


20dBm 


Noise power spectral density 


-174dBm/Hz 



A. Channel Model 

The simulations are implemented using the COST-231 channel model as described in [15]. The model 
assumes both the BS and relays are at some height off the ground and treats the BS-relay channel as 
Rician. The BS-destination and relay-destination channels are modeled as Rayleigh. The path loss in the 
BS-relay channel is made up of two components, free space loss and multi-screen loss. In addition to these 
two, the BS-user and the relay-user channels have a roof top-to- street diffraction loss. For the values of the 
parameters that we consider, the COST-23 1 channel model suggets a distance attenuation in channel power 
of 20dB/km for the first 657 meters and 38dB/km for greater distances. The model therefore appears to 
be conservative in the sense that one would expect the LoS component in the Rician fading to attenuate 
slower than the other non-LoS components. In the MISO case, the large scale fading in all the channels 
between the transmit antennas and a particular user, is the same. Each user is assigned an orthogonal 
channel of bandwidth of 200kHz, resulting in a noise power of -120dBm. The chosen system parameters 
are given in Table HI 

B. Decoding at the relays 

To form a tractable problem, we had made the assumption that the relays always successfully decode 
the message transmitted by the BS in the first time slot, and the data rate is the limited by compound 
source-relay-destination channel capacity, as in Eqn. Q. To verify the assumption, we consider a circular 
cell, centered at a BS, of radius one kilometer with J = 4 relays positioned at (±200\/2m, ±200v^m), 
i.e., on a ring of radius 400m. 3x10^ user locations in the cell are randomly generated. For each location, 
independent channels are generated using the channel model. As shown in Fig. [H we divide the cell into 
annular rings of radius 100 meters. In Table |ll] we list the percentage of number of locations where 
Eqn. Q is valid. It is evident from the table that the assumption we make is valid whenever the user 
is farther than 300m from the BS. Essentially, for all user locations of interest, i.e., areas where users 
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TABLE II 

Percentage of locations where Eqn. i|6j is satisfied 



Distance from 


% 


Distance from 


% 


the BS (m) 


locations 


the BS (m) 


locations 


0-100 


93.591 


500-600 


99.943 


100-200 


99.642 


600-700 


99.963 


200-300 


99.815 


700-800 


99.977 


300-400 


99.309 


800-900 


99.989 


400-500 


99.482 


900-1000 


99.992 



have a relatively weak channel to the BS, the assumption is valid. It is worth emphasizing that these are 
conservative numbers. 

C. Tightness of the bounds 

Our next of simulations test the tightness of the upper bound as developed in this paper and the resulting 
heuristic which acts as a lower bound. Note that this heuristic is our final solution to the joint selection 
and power allocation problem. The relay assignment and the power allocation is done based on the 
instantaneous channel powers. For this simulation, the channels are generated as independent realizations 
of a unit-variance Rayleigh fading random variable. For a fair comparison, the power allocated to each 
relay is set to 1/J, i.e., all curves use the same total power. The curves presented here are averages over 
one thousand random user locations. 

Figure |2] plots the upper bound, and the sum rate achievable by the heuristic (that also acts as a lower 
bound on the achievable sum rate) for varying values of J and K. The average total transmit power to noise 
power ratio is set to 30dB. As is clear from the figure, the upper and lower bounds are indistinguishable. 
As explained in Section Hill this is because it is quite rare for a user to be allocated power from multiple 
relays, i.e., selection is essentially inherent in the approximate solution. The heuristic, therefore, is an 
extremely effective solution to the joint selection and power allocation problem. By an exhaustive search, 
we also find the exact maximum sum rate for the case with J = 2 relays and K between one and eight. 
Note that since each exhaustive search requires solution of water-filling problems, any larger value of 
J is infeasible. 

Figure [3] plots the upper and lower bound to the max-min rate for varying values of J and K averaged 
over many channel realizations. In this simulation, the average total transmit power to noise power ratio 
is set to 20dB. Again, the bounds are extremely tight and the heuristic provides an effective solution. The 
slight difference is due to the rare case where a user is allocated power by two relays (see Section HIH) . 
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In interests of brevity, we do not provide a similar plot for the max-sum rate with a rate constraint. 

D. Results for a cellular network 

In this section, we use the theory developed for solving the max-min and max sum rate problems, to 
estimate the performance gains, with respect to a SISO and MISO system in cellular network setting. We 
consider a cell of radius rcdi. Performance here is measured as the increase in cell-size made possible by 
relaying. Since we wish to study the improvement in the rates to the users with poor channels to the BS, 
we consider users in the outer annular ring, of inner radius rccii/2 and outer radius iceii, the area shaded 
in gray in Figured] Users are distributed uniformly in the region with a constant user density of (SO/vr) 
per square kilometer. We consider the following three system models for comparison: 

1) A cellular network with a single antenna BS, communicating to multiple users with single antenna 
receivers (multiuser SISO system). 

2) A cellular network with a BS with five transmit antennas, communicating to multiple users with 
single antenna receivers (multiuser MISO system). 

3) A cellular network with a BS with a single antenna and assisted by four relays positioned on a ring 
of radius 0.4rceii, communicating to multiple users with single antenna receivers. 

For the simulation, we generate 50 random sets of locations for the users. We then use the COST-231 
model to generate the BS-user and relay-user channels. For each set of locations, we generate one set 
of large-scale fading variables. To average over small-scale fading random variables, for every set of 
locations, we generate 500 small-scale fading random variables. As indicated in Table HI the total power 
used in communication is set to 20dBm. 

In the first example, powers are allocated to maximize the sum rate to all the users. For a fair comparison, 
we use this to compute the data rate averaged over all users. In the SISO case, the system uses water-filling 
to allocate power to the multiple users. In the MISO case, the BS is assumed to know the channel vector 
to each user and can both match to the channel and allocate power via water-filling. Finally, in the case 
with relays, selection and power allocation uses the scheme developed in Section [nil 

In Figure IH we plot the average user rate as a function of the radius of the cell. We compute the rates 
as given by the lower bound, assuming that the power allocation is done only in the second time slot. 
In the first time slot, the BS distributes the power equally among all the users. This is done to ensure 
that the relays are able to decode all the transmitted messages. In the second time slot, each of the relays 
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uses one fourth of the available power to communicate with the users it assists. This ensures that the 
total power used is the same in all the three system models. Interestingly, Figure |4] shows that a MISO 
system provides higher average data rates (and hence the sum rate) compared to the system with relays. 
We explain these graphs in the following section. 

Next, we repeat the simulations by allocating power using the max-min algorithm, and then computing 
the outage rates for each of the system models. For the SISO and MISO cases, computing this power 
allocation is fairly straightforward. The optimal power allocation is the one such that all the users have 
the same data rate. When relays are employed, we use the methodology developed in this paper to solve 
the max-min optimization problem. 

We plot the outage rates for 10% and 1% outage, as a function of rcdi in Figure [51 Here we see 
a reversal in performances, with the system with relays providing higher outage rates compared to the 
MISO system. As expected, the BS-user communication in these systems is more susceptible to channel 
fluctuations. This plot is discussed further below. 

E. Discussion 

A user in a heavily shadowed region has a weak channel to the base station. Having multiple antennas 
at the base station does not help much. Relays aid such users by providing alternate paths to the base 
station. This is consistent with the data in Figure [5] where a relay system provides higher outage rates. 
This is because the outage rates depend on the data rates to the users with weakest channels. On the other 
hand, the max-sum rate algorithm, allocates more power to the users with the strongest channels. The 
MISO system provides higher data rates compared to a relay system. As shown in Fig. H the loss in half 
the bandwidth incurred in switching from direct transmission to co-operative transmission outweighs the 
benefits brought by the additional diversity. 

In a network setting where a user has the same average channel to all the J relays, selection cooperation 
achieves order J + 1 diversity [6]. However, because of the geometry of a cellular network and because 
of the rapid deterioration of the channel powers with distance, most users have good channels to only a 
small set of relays. The ejfective diversity order is, therefore, limited. 

Figures H] and [5] also lead to a cautionary result. These results indicate that, compared to a system with 
SISO communication, deploying relays does offer substantial improvements. The area serviced effectively 
by a single BS, helped by relays, can significantly expand. However, these improvements need to be 
commensurate with the infrastructure costs involved in the deployment of these relays including both the 
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antenna system cost and 'non-technical' costs such as the required real-estate. If the cost of a relay were 
on the same order of magnitude as a base station, the improvements in the cell radius, as shown by the 
simulations do not justify the additional cost. Also, depending on the performance metric, a MISO system 
may perform better or almost as good as the relay system, but with significantly lower costs. 

Clearly, a complete financial cost/benefit analysis is beyond the scope of this paper. Furthermore, the 
examples presented here are limited and do not explore every potential parameter. However, do note that 
the our results are optimistic in assuming the relays can always decode and that the transmitters have all 
the information they need to make optimal decisions. Our goal here is to indicate that significant gains are 
possible, but are context and scenario dependent. These results also indicate the need for exploring alternate 
ways of exploiting cooperative diversity. We also need to explore alternate hybrid schemes wherein the 
relays help only those users who need it. 

V. Conclusion 

This paper deals with the use of cooperation in a cellular network wherein a base station is assisted by a 
few dedicated relays. Previous work largely for mesh networks has shown the importance of selection, i.e., 
each user using only one relay, since this minimizes the overhead due to orthogonal channels. However, 
in a scenario with multiple data flows, selection has been either brute force or ad-hoc. Previous work 
has also largely ignored the problem of power allocation once the selection is achieved. In this paper we 
developed an optimization framework to solve the problem of joint selection and power allocation. 

The optimization problem uses the achievable sum rate and max-min user-rate as two figures of merit. 
Given that the selection problem has exponential complexity, in this paper we formulate alternative convex 
optimization problems whose solution provides upper bounds on the two metrics. However, for practical 
values of number of users, the bound is indistinguishable from the true solution. Since this solution can 
violate the selection condition, a related heuristic is derived that assigns users to the relay which allocates 
it the maximum power. The resulting lower bound is also extremely tight and we have an efficient solution 
to the problem at hand. The numerical examples, using realistic channel models, illustrate the benefits 
achievable due to relaying. 
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Fig. 1. A relay aided cellular network 



An upper bound(UB) and a lower bound(LB) to the maximum sum rate 
for different number of relays(J) 




5 10 15 20 25 30 35 40 

Number of users(K) 



Fig. 2. The proposed upper bound to the maximum sum rate and the heuristic (a lower bound) as a function of the number of users, 
that both the bounds are extremely tight. 



Note 
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Fig. 3. The proposed upper and the lowerbound to the max-min rate. 




Fig. 4. 



Average data rate as a function of the radius of the cell. 
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Fig. 5. Outage rate as a function of tiie radius of the cell. 



